o 

(N 

> 

O 



cr 



(N 
> 

O 

ON 
(N 

l> 

O 



X 



Extending the effective-one-body Hamiltonian of black-hole binaries to include 
next-to-next-to-leading spin-orbit couplings 
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In the effective-one-body (EOB) approach the dynamics of two compact objects of masses mi and 
?ri2 and spins Si and S2 is mapped into the dynamics of one test particle of mass /i = mi 7712/(7711 + 
TO2) and spin S, moving in a deformed Kerr metric with mass M = mi + 777,2 and spin Skcit- In 
a previous paper we computed an EOB Hamiltonian for spinning black-hole binaries that (i) when 
expanded in post-Newtonian orders, reproduces the leading order spin-spin coupling and the leading 
and next-to-leading order spin-orbit couplings for any mass ratio, and (iii) reproduces all spin-orbit 
couplings in the test-particle limit. Here we extend this EOB Hamiltonian to include next-to- next- 
to-leading spin-orbit couplings for any mass ratio. We discuss two classes of EOB Hamiltonians 
that differ by the way the spin variables are mapped between the effective and real descriptions. 
We also investigate the main features of the dynamics when the motion is equatorial, such as the 
existence of the innermost stable circular orbit and of a peak in the orbital frequency during the 
plunge subsequent to the inspiral. 

PACS numbers: 04.25. D-, 04.25.dg, 04.25.Nx, 04.30.-w 



I. INTRODUCTION 

Coalescing compact binaries composed of neutron 
stars and/or black holes are among the most promis- 
ing gravitational-wave sources for ground-based detec- 
tors, such as the Laser Interferometer Gravitational- wave 
Observatory (LIGO) [1], Virgo [2], GEO [3], the Large 
Cryogenic Gravitational Telescope (LCGT) [4], and fu- 
ture space-based detectors. 

So far, the search for gravitational waves with LIGO, 
GEO and Virgo detectors has focused on non-spinning 
compact binaries [5-9], although in Ref. [10] single-spin 
templates were used to search for inspiraling spinning 
compact objects. Within the next 4-5 years LIGO and 
Virgo detectors will be upgraded to a sensitivity such 
that event rates for coalescing binary systems will in- 
crease by a factor of one thousand. Thus, it is timely 
and necessary to develop more accurate templates that 
include spin effects. For maximally spinning objects, we 
expect that reasonably accurate templates would need 
to be computed at least through 3.5PN order. In the 
non-spinning case, studies at the interface between nu- 
merical and analytical relativity have demonstrated that 
templates computed at 3.5PN order are indeed reason- 
ably accurate. 

In the last few years, motivated by the search for grav- 
itational waves, the knowledge of spin effects in the two- 
body dynamics and gravitational-wave emission within 
the post-Newtonian (PN) 1 approximation has improved 
considerably. In particular, spin-orbit (SO) effects in 
the two-body equations of motion are currently known 



1 We refer to nPN as the order equivalent to terms 0(c 2 ") in the 
equations of motion beyond the Newtonian acceleration. 



through 3.5PN order (i.e., 2PN order beyond the lead- 
ing SO term) [11-19], and in the energy flux through 
3PN order [12, 13, 20-23] (i.e., 1.5PN order beyond the 
leading SO term). Moreover, spin-spin (SS) effects have 
been computed through 3PN order (i.e., 1PN order be- 
yond the leading SS term) in the conservative dynam- 
ics [11, 20, 24-34] and also in the multipole moments [35]. 

In order to build reliable templates and search for 
gravitational-waves from high-mass compact binaries 
that merge in the detector bandwidth, it is crucial to 
improve the PN approximation by resumming the dy- 
namics and gravitational emission in a suitable way and 
by using numerical relativity and perturbation theory as 
a guidance. The effective-one-body approach (EOB) [36- 
40] offers the possibility of fulfilling this goal. The EOB 
approach uses the results of PN theory, not in their orig- 
inal Taylor-expanded form (i.e., as polynomials in v/c), 
but instead in a suitably resummed form. In particular, 
it maps the dynamics of two compact objects of masses 
777i and 7772, and spins Si and S2, into the dynamics 
of one test particle of mass \x = m\ 1712/(1711 + m 2 ) and 
spin S* moving in a deformed Kerr metric with mass 
M = mi +7712 and spin Sscerr- The deformation parame- 
ter is the symmetric mass ratio 77 = 7771 7772/(7771 + TO2) 2 , 
which ranges between (test particle limit) and 1/4 
(equal- mass limit). The analyses and theoretical progress 
made in Refs. [41-57] have demonstrated that faithful 
EOB templates describing the full signal (i.e., the inspi- 
ral, merger and ringdown) can be built and used in real 
searches [9]. 

Here we build on previous work [39, 46, 58, 59], employ 
the recent results of Ref. [19] and extend the EOB con- 
servative dynamics, i.e. the EOB Hamiltonian, through 
3.5PN order in the SO couplings. Since the mapping 
between the PN-cxpanded Hamiltonian (or real Hamil- 
tonian) and the EOB Hamiltonian is not unique, we ex- 



plore two specific classes of EOB Hamiltonians, which 
differ by the way the spin variables of the real and effec- 
tive descriptions are mapped. 

This paper is organized as follows. In Sec. II, af- 
ter reviewing the logic underpinning the construction 
of the EOB Hamiltonian, we proceed in steps and ex- 
tend the EOB Hamiltonian proposed in Ref. [46] through 
3.5PN order in the SO couplings. In particular, in 
Sec. II A we derive the PN-expandcd Arnowitt-Dcscr- 
Misncr (ADM) Hamiltonian in the EOB canonical co- 
ordinates; then, after computing in Sec. II B the effec- 
tive Hamiltonian corresponding to the canonically trans- 
formed PN-expanded ADM Hamiltonian, we compare it 
(in Sec. II C) to the deformcd-Kcrr Hamiltonian for a 
spinning test-particle [46], and work out (in Sees. II D 
and HE) two classes of EOB Hamiltonians. In Sec. Ill 
we study the dynamics of these Hamiltonians for equa- 
torial orbits, and in Sec. IV we summarize our main con- 
clusions. 

We use geometric units G = c = 1 throughout the 
paper, except when performing PN expansions, where 
powers of the speed of light c are restored and play the 
role of book-keeping parameters. 



II. THE EFFECTIVE-ONE-BODY 

HAMILTONIAN FOR TWO SPINNING BLACK 

HOLES 

The main ingredient of the EOB approach is the 
real PN-expandcd ADM Hamiltonian (or real Hamilto- 
nian) describing two black holes with masses 7711,7712 
and spins Si, S2. The real Hamiltonian is then canoni- 
cally transformed and subsequently mapped to an effec- 
tive Hamiltonian H c g describing a test-particle of mass 
/j = mi 7772/ (mi + 7772) and suitable spin S*, moving in 
a deformed Kerr metric of mass M = m\ + 7772 and suit- 
able spin Skcii- The deformation is regulated by the 
binary's symmetric mass-ratio parameter, 77 = /jl/M, and 
therefore disappears in the test-particle limit 77 — ^ 0. The 
so-called improved real (or EOB) Hamiltonian reads 



rri in proved 
-"real 




(1) 



The computation of the EOB Hamiltonian consists of 
several steps. We briefly review these steps and the un- 
derlying logic that we will follow in the next sections: 

(i) We apply a canonical transformation to the PN- 
cxpanded ADM Hamiltonian using the Lie method, 
obtaining the PN-expandcd Hamiltonian in EOB 
canonical coordinates (see Sec. II A); 

(ii) We compute the effective Hamiltonian correspond- 
ing to the canonically transformed PN-expandcd 
ADM Hamiltonian (see Sec. II B); 



(iii) We PN-expand the deformcd-Kcrr Hamiltonian for 
a spinning test-particle derived in Ref. [46] (see 
Sec. II C); 

(iv) We compare (ii) and (iii), and work out the map- 
ping between the spin variables in the real and 
effective descriptions, and compute the improved 
EOB Hamiltonian (see Sees. II D and HE). 



A. The ADM Hamiltonian canonically transformed 
to EOB coordinates 

Following Ref. [46] , we denote the ADM canonical vari- 
ables in the binary's center-of-mass frame with r' and p' , 
and we introduce the following spin variables: 



<T = 


= S!+S 2 , 


(2) 
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Henceforth, to keep track of the PN orders, we rescale 
the spins variables as cr* — >• cr* c and cr — > cr c. 

We use the spin-independent part of the ADM Hamil- 
tonian through 3PN order [38], and we include SO ef- 
fects through 2PN order beyond the leading-order effects 
(1.5PN), thus through 3.5PN order. In particular, the 
ADM SO Hamiltonian at 3.5PN order was computed re- 
cently in Ref. [19] (the ADM SO Hamiltonian at 1.5PN 
was computed in Ref. [60], and at 2.5PN in Ref. [58]). In 
the binary's center-of-mass, the ADM SO Hamiltonian 
reads 
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where we indicate L' = r' x p' and 



ADM 



= 2 + - 
Cr 



19 -, 

-7TVP 



^(n'-p'Y 



2r i)-r 



^fin'-P'f 



21 



(I + 77) 



+ \r, (-9 + 2277) p' 4 - -U (314 + 3977)^ 
8 16 r 



P' 2 



1 
16 



M 



T) (256 + 45t;) — (ri ■ p') 



'\2 



.± v (-A + 9 V )p' 2 {n'-p') 2 



(5a) 



ADM 
9a* 



2r ] )p' 2 



-T) (n' ■ p') 2 



-(5 + 277)^ 



1(75 + 8277) (M. 



--(7-37r ? + 397y 2 )p' 4 

16 

~ V (32 + 15r,)^-(n'-p, 
16 r 

±r,{-l + 2r,)p l2 {n'-p'f 



r 



(5b) 



with n' = r'/r', and where we have introduced the 
rescalcd conjugate momentum p' = p' / fi. 

In order to canonically transform the ADM Hamilto- 
nian to EOB coordinates, various approaches are possi- 
ble. A popular method, used in the previous work on the 
EOB model [37, 38, 58], is to use a generating function 
that produces a near-identity transformation, i.e. one of 
the form G(q', tt) — q n iTi + e G(q', tt), where (q, it) are the 
phase variables (including the angles defining the spins 
and their conjugate momenta, see Ref. [59]) and e is a 
small parameter. Expressing the initial "primed" coor- 
dinates (the ADM coordinates) in terms of the new "un- 
primed" coordinates (the EOB coordinates), one gets 



dG(q',7T) , 

- e — - = q - c - 

OTTi 

c2 d 2 G(q,7T) dG(q,7T) 
dmdqi diTj 
dG(q',TT) 
+ e ^q^ =7Tl 
c2 d 2 G(q,7T) dG(g,ir) 
dq i dqi dnj 



dG(q,7T) 
dir l 

0(e 3 ), 

dG(q,iT) 
dq i 

0(e 3 ). 



(6a) 



(6b) 



Because under a time-independent canonical transforma- 
tion the Hamiltonian transforms as H(q,p) = H'(q',p'), 
Eqs. (6) imply 



H(q,Tr) ---- H'(q',7T')=H'(q,7T) 
~dH'(q,iT) dG{q,ir) 



+r 
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dq l diTj 
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+ 0(e 3 ). (7) 



The terms of order 0(e) in this equation can be rewritten 
as e{G, H'}, which is very convenient because it trans- 
forms a sum over all the phase variables (including the 
angles defining the spins and their conjugate momenta) 
into a Poisson bracket that can be computed using only 



the commutation relations {x l ,pj} = <K, {x l , S 3 , a s} = 0, 

{Pi> S la)} = °> and ^(o)'^)} = 5 (*){b)^kS k (a) {a,b = 
1 , 2 being indices that distinguish between the two black 
holes). Unfortunately, the terms 0(e 2 ) cannot be eas- 
ily expressed in terms of Poisson brackets, which makes 
them hard to compute (because the spin variables must 
be carefully taken into account in the sums). Also, the 
generalization of Eq. (7) to higher orders in e becomes 
more and more complicated. 

A possible alternative to the generating function 
method mentioned above is given by the so-called Lie 
method [61]. This approach exploits the fact that the 
flux of the Hamilton equations is canonical. Therefore, 
one can define a fictitious Hamiltonian T-L{q, it) whose flux 
sends some initial data {q,ir) to (q'(q,m e),ir'(q, tt, e)), 
where e is the "time" variable of this fictitious Hamilto- 
nian. The canonical transformation is then simply given 
by (q'(q, tt, e), Tr'(q, m e)). The advantage of this approach 
is that any function f(q,ir) satisfies / = {/, 1-L} (where 
we denote with ' = d/de) under the Hamiltonian flux of 
W. Defining for convenience Q = —%, this equation be- 
comes / = {Q, /}, and denoting the differential operator 
{Q, . . .} by Lq, a Taylor expansion yields 

oo n 

f(q'(q,me),7T'(q,7T,e)) = V - £g f{q,rr) 

= exp(eCg)f(q,Tr) 

= f(q,7T)+e{g,f}(q,7T) 



1 



£ 2 {S,{S,/}}(^) + 0(£ 3 ). 

(8) 



Specializing to the (non-fictitious) Hamiltonian H' , we 
obtain the equivalent of Eq. (7), that is 

H(q,7T) = H'(q'y)=H'(q,iT)+e{g,H'}(q,iT) 

+ \e 2 {g,{g,H'}}{q,ir) + 0{e*). (9) 

As already mentioned, the above expression allows us to 
account for the spin variables very easily, if necessary 2 , 
by means of the commutation relations {x l ,Sj a A = 0, 

{Pi> S (a)} = °i and { S U> S (b)} = S (o,)(b)eijkSf a) . 

In this paper we will use the Lie method to generate 
the canonical transformation from ADM to EOB coordi- 
nates. In particular, we assume 

G(r,p) = r-p + g NS (r,p) + g s (r,p,(T*,cT), (10) 

where (?ns is the purely orbital part of the fictitious 
Hamiltonian, while t/s is the spin-dependent part, which 



2 The Poisson brackets of the spin variables with themselves do 
not enter in the computations that we perform in this paper, but 
they do enter at higher PN orders. 



we assume to be linear in the spins since in this paper we 
focus on the SO terms only. Because the transformations 
(7) and (9) agree at leading order in the perturbative pa- 
rameter e, G and 5 must agree at leading PN order. In 
particular, since the purely orbital generating function 
for the transformation from ADM to EOB coordinates 
starts at 1PN, 5ns must start at 1PN order too, that is 

5ns (r,p) = 5nsipn(7 , ,p)+5ns2Pn(t , ,p)+0 ( -g J , 

(11) 

where 5ns ipn must coincide with Gns ipn , and therefore 
be given by [37] 



1 



5nsipn(t",p) = -j(r-p) 



1 2 M f 1 
._„p +-(l + -r/ 



(12) 

At 2PN, instead, 5ns does not coincide with Gns, but a 
computation similar to the one in Ref. [37] easily shows 
that 



5nS2Pn(»",p) 
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,(13) 



with 



a = — , 



= {(*-ri), (14a) 



4 + 77 1-777 + y 2 
7 = V—z—, 6= . (14b) 



[Note that the functional form (13) is the same as for 
Gns2PN, but the values of the parameters a, j3, 7 and <5 
are different from those of Ref. [37].] 

Similarly, the spin-dependent part of the fictitious 
Hamiltonian, 5s , must start like G$ at 2.5PN order: 



5s(r-,p,cr*,cr) = 5 S 2.5PN(r,:P,er*,cr) 



+ 5s3.5PN(r,p, C r*,<T)+0( -) . (15) 



and if we restrict to functions that are linear in the spin 
variables, it must be [15] 



5s2.5PN(r,p,<T*,cr) = 
1 



(r-p) [a {r,){L-(T) + b {ri){L-<T*)] , (16) 
c w 1 ~ 

where 00(77) and &o( 7 7) are arbitrary gauge functions. 
[Note that restricting to functions that are linear in the 
spin variables is justified because here we are looking at 
SO effects only, but in general cubic terms in the spin 
may be present, see Ref. [46].] 



The most general form for 5s at 3.5PN order is instead, 
if we restrict again to functions linear in the spins, 

5s3.5PN(r,p) = -^-3 (r • p) {(L ■ er) [01(77) p 2 

, M ... ,2 

+02(77) \-a 3 (r}) (n-p) 

r 



+ {L.a*) 



bi(il)p 2 +b 2 (r]) — 



i(r])(n-p) 2 J 



(17) 



where ^(77) and bi(rf) with i = 1,2,3 are other arbitrary 
gauge functions. To ease the notation, henceforth we 
drop the 77 dependence in the gauge parameters, both 
at 2.5PN and 3.5PN order, and will denote them simply 
with a, and bi (with i = 0,3). 

Applying Eq. (9), we obtain that the 3.5 SO Hamilto- 
nian in EOB coordinates is given by 

tt _ rrADM , in rrADM l 

«S03.5PN - -«S0 3.5PN + I^NS2PN,^soi.5PN} 

+ {£nS1PN,# S 2.5Pn} 

1 If. rrADMl 

+ tt/S2.5PN5-n 1PN i 
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(18) 



A tedious but straightforward calculation gives the sev- 
eral terms entering the above equation: 
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Also, we have [46] 



Rcf. [46], while 
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where #£™, E ivT can bc found in Rcf - [37], 
the explicit expressions of {£s2 spn, -^ncw?} and 
{0ipn,H£?}W arc given in Eqs. (5.20), (5.24) of 
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Also, we note that Eqs. (26b) and (26d) immediately im- 
ply S2.5PN = G2.5PN, i-e. the 2.5 PN gauge parameters 
ao and &o appearing in Eq. (16) have the same mean- 
ing in the Lie method and in the generating function 
approaches. 



B. The spin-orbit terms in the effective 
Hamiltonian through 3.5PN order 
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Following Refs. [37, 39, 62], we map the effective and 
real two-body Hamiltonians as 
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where -ff rea i is the real two-body Hamiltonian contain- 
ing also the rest-mass contribution Mc 2 . Expanding 
Eq. (28) in powers of 1/c, we have 



Using Eqs. (18) and (26), we find that through 3.5PN 
order the SO couplings of the effective Hamiltonian are 
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C. The PN-expanded Hamiltonian of a spinning 
test-particle in a deformed Kerr spacetime 
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The deformed-Kerr metric was obtained in Ref. [46], 
and it reads 
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The potentials in these equations are given by 
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n 9 1 M 3 
cDfd = 2a M r + a 77 w fd M + a 77 a; fd , 



(33) 

(34) 
(35) 
(36) 



(37) 



where wj? d and Wf d are two "frame-dragging" parameters 
(that we will fix later), and where 



2M 2nM 3 /94 41 9 \ n 



D^(r) = 1 + -^ + 2(26-3,?) ^ r - 



Af 4 



(38a) 
(38b) 



The Hamiltonian of a spinning test-particle in the 
dcformcd-Kcrr spacetime is 



and 



ff = 2- 



ij 9 U 9 tJ 



Hq — H : 



so 



H, 



ss 



7 

(42) 
(43) 



(44) 



H = H NS + H S , 



(39) 



with 



#NS =/3 i p i + «i/m 2 + 7 ij >iPj + Qa(p) , (40) 

where the term Q4 (p) is quartic in the space momenta pi 
and was introduced in Ref. [62]. Moreover, we have 



a 



V^' 



(41) 



where -ffso contains the odd terms in the spins (and 
therefore, in particular, the SO terms) and -ffgg contains 
the even terms in the spins (and therefore, in particular, 
the spin-spin terms of the kind SkottS*). 

Since here we are interested in the SO couplings, we 
consider only Hso (Hss can be read from Eq. (4.19) in 
Ref. [46]): 



Hso = 



2u ~ fi (eP+'-B) (p-£r)(S-S Ko 



="-2£ 



-fJ-cos e (p ■ n) £ 2 - ^/Q (y r (p ■ v r) + (^ cos e - ^cos e) (P • «) £ 2 ) 



(S-£)J » r (p-vr)(y/Q + r 
(P-Sr) (2y/Q + l) 



B 2 + d"+ v 



Jv r (S-v)-v cos0 (S-n)^ B-JBre^lp-tr) [VQ + 1 (5 • v) 



(45) 



r 



Where S K err = SKerr/Skcrr, I = ^Kerr X 71, V = 71 X £, 

and where 



1 + 



A r (p-n) 2 (p-|r) 2 £ {p-vrf 



A t sin 2 i 



E sin 2 9 



and 



r (r 2 + a 2 ) \(r 2 + a 2 ) A' t - 4r A t l 

u r = - + 



V CO s8 



Mr- 



2A t A t 

a 2 (r 2 + a 2 ) cos 6 (r 2 + a 2 - A t ) 



/AT' 



A*E 



(J-COS0 



a 2 cos 8 



B 



B r = 



/AZAi -2 A* 



2fi _ y 2u _ A * S 

e _S, e - 



2 V / A7A7 

, J = ^/aT, 



(46) 

(47a) 
(47b) 
(47c) 
(47d) 
(47e) 



in which we use a prime to denote the derivative with 
respect to r. To obtain the SO couplings through 3.5PN 



order, we expand Eq. (39). In particular, it is suffi- 
cient to consider the first term in the right-hand-side of 
Eq. (40), and set a = (dcformcd-Schwarzschild limit) in 
Eqs. (45), (46) and (47). Doing so, for the PN-expanded 
deformed-Kerr Hamiltonian we obtain 

' L ■ OKcrr j 



H, 



NS 

S0 1.5PN 



H, 



NS 

S0 2.5PN 



r 3 c 3 

— jjWfd — L-S Ke 



rrNS 
-"SOS.SPN 



1 



■»?w fd 



M 



L S 



Kerr • 



and 



H, 



1.5PN - ^J 

1 



L S* , 



(48a) 
(48b) 

(48c) 
(49a) 



-"302. 5PN 



H, 



'3.5PN - ^7 



L-S* 



L-S* 
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(49b) 



i(-l-42 ?/ + 6r, 2 ) (¥- 



7 ~4 1 + 

— v -\ 



G// 



V 



{n-pf 



(49c) 



D. The EOB Hamiltonian: spin-mapping 
dependent on dynamical variables 

We now determine the mapping between the spins cr 
and cr* of the effective ADM Hamiltonian and the spins 
■Skcit and S* of the EOB Hamiltonian by imposing that 
the deformed-Kerr Hamiltonian given by Eqs. (48) and 
(49) coincides with the effective Hamiltonian given by 
Eqs. (30) and (31). As found in Ref. [46], we have to 
assume that the mapping depends on the orbital dynam- 
ical variables p 2 , n • p and r. The general mapping of the 
spins has the form 



S* = 



- A« + - A (2 J 
c 2 a c 4 " 



(50a) 



Skc- = rr + 1 A « + 1 A« 



(50b) 



At 2.5PN order, if we assume CJf d = [see Eq. (37)] and 
A ff 1} = 0, we have [46] 



(i) 



i(-4fe + 7r ? )^ + i(26 + ?? )(g-l) 



2 Ml 

-^(ao + r/) — + -(8a + 3r y )(g-l) 

3 r 12 



-(2a + 37 1 )^(n-p) 2 



(51) 



and at 3.5PN order, assuming Wjy = [see Eq. (37)] and 

(2) 

Act ; = 0, we obtain 
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A\2 



1 



(— 14ao — 6a 2 — 56t/ — 15ao?y — 21?/ 2 ) 



^(-16o 3 + 24a 7? + 27 ?? 2 ) |f (n • p) 4 



144 



H 



-96ai - 45r/ - 24a 7/) (Q - l) 2 + — (-14a - 24 fll + 24a 2 - 10971 + 66a 7/ + 51t/ 2 ) — (Q - 1) + — (4a 

36 r 24 



-48ai + 48a 3 - 
-1477/ 2 )Ml| 



6r/ - 60aov - 39 ? / 2 ) ^ (n • p) 2 (Q - 1) + 1 (- 



(n-p) ; 



32ai - 64a 2 - 48a 3 - 16?/ - 48a ?/ 



(52) 



Note that as in Ref. [46] , we have replaced, in the expres- 
sions for A CT „ and A CT « , the term p 2 with j^Pipj = Q — 1 
and the term (n ■ p) 2 with A r (n ■ p) 2 /S = </ rr /3 2 . 

Having determined the spin mappings, wc can write 
down the real improved (or EOB) Hamiltonian for spin- 
ning black holes, which turns out to be 



rrimprovcd 
-"real 



MW 1 + 277 



H, 



cir 



(53) 



where 



H s +f3 l Pi + a Jfj? + rfi Pi p 3 + Q4 (p) 



/'• 



2Mr 



(<**■ 



3n*ra J )S;s; 



(54) 



E. The EOB Hamiltonian: spin-mapping 
independent of dynamical variables 

In the previous section, we had to assume a dependence 
on the orbital dynamical variables p 2 , n • p and r in the 
mapping between the spins cr and cr* of the effective 
ADM Hamiltonian and the spins Sjcerr anc ^ &* of the 
dcformcd-Kcrr Hamiltonian. To avoid this dependence 



on the dynamical variables and obtain the much simpler 
mapping 



S* = a* 

^Kerr = O" ■ 



(55a) 
(55b) 



we need to modify the Hamilton- Jacobi equation by 
adding terms depending on the momenta and spins. 
Since in this paper we are dealing only with SO effects, we 
will neglect modifications that involve spin-spin terms. 
We start by observing that in the presence of spins the 
linear momentum P M , which is related to the canoni- 
cal momentum by P p = p p + EP a S* pa [see Eq. (3.28) 
of Ref. [59]], satisfies the Hamilton- Jacobi equation [59] 



l^+P^P" 



fi +P^P"+ 



2^/s;- 



■0(S* 



0. (56) 



Here, S* is the spin tensor of the test-particle [see 
Ref. [59] and also Eqs. (2.4)-(2.7) in Ref. [14]]. Equa- 
tion (56) leads to the correct Hamiltonian for a spinning 
particle in curved spacetime, at linear order in the parti- 
cle's spin [59]. To modify the Hamilton- Jacobi equation, 
a suitable ansatz is 

/i 2 + 9 ^(V, S K err)^ P„ + 2E^pf i S* pa 

+ [B^ x (r)p pPl/ p x 

+ B% XTa (r)p,p„p x p T p a ]Sr 

+ A" uXT (r,S Ke „)p^PuPxPr 

+ A^ XTpa (r, S Km )p^PuPXPrPpPa 

+ • • • = . (57) 

If we make use at lowest order of the condition /i 2 + 
g^fiPp,Pv — 0, we can replace pt with the spatial compo- 
nents of the momentum, and obtain the following gener- 
alized form of the effective Hamiltonian 

H c g = 

/3 l Pl +a sjy 2 + yipiPj + Qa{p) + Qs{r,p, S*,S Ke ,r) 



/' 



2Mr 



(5 ij - 3n* n j ) S* S^ + H s 



(58) 



where Qa{p) is a quartic term in the momenta [62], 
which is due to the presence of the quartic term 
A^ XT Ptl p v pxPr in Eq. (57), and Qs(r,p,S*,S KeiI ) is 
a term linear in S* and Skcit 



■)5k. 



S ' i 



Qs(r,p,S*,S KeII ) = Q? K °"(r,p)S* Ke „ + Qf (r,p)Sl . 

( 59 ) 
In particular, the term Qf {r,p)Sl comes from the 

terms B^ x {r) Ptl p,p x S pa and B% x ™ p p p v p x p r p a S pG 

in Eq. (57), while the term Qf Ko " (? , ,p)5'kcit comGS from 
A tyXr PtPuPxPr and A tvXTpa p t p v p x p T p p p a (through 
the dependence of the tensors A pvXr and A pvXrpcr on 
■Skerr)- Finally, the term Hg in Eq. (58) comes from the 
presence of 2E p,7 p p S* in Eq. (57). As already stressed, 
this happens because Eq. (56) leads to the correct Hamil- 
tonian for a spinning particle in curved spacetime, and 
in particular to Hs, which is the spin-dependent part of 
that Hamiltonian [59]. 

Through 3.5PN order the quantities Qf Ko "(r,p) S^ crl . 
and Qf (r, p) SI must have the form 



QKr,py 



A* i k i 

273 £ ^ n > S X 



r c 
1 / M 



-o I Ci h c 2 p + c 3 (n • p)' 



C4P 4 + c 5 [ — ) 4- c 6 (n • ft) 4 



c 7 p — + c 8 (n • pY — + c 9 (n • p p z 
r r 



(60) 



where s stands for either Skcit or S* , while the coeffi- 
cients c n , n = 1 , ... 9 are determined by the mapping of 
the effective to the real description. A straightforward 
computation leads to 



Qs = Qs2.5PN + Qs 3. 5PN 



(61) 



where 



QS 2.5PN 0" , P, S* , S K err) 



(S'Kerr - L) 



-(S*-L) 
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I r I Zj 



. (62) 



Qs3.5PN(f, p, S*, SKerr ; 
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4 r 2j 



(63) 



r 



Finally, the EOB Hamiltonian is obtained by inserting 
Eq. (58) into Eq. (53). 



III. THE EFFECTIVE-ONE-BODY DYNAMICS 
FOR EQUATORIAL ORBITS 

We stress that the EOB models introduced in the pre- 
vious sections have the correct test-particle limit, for 
both non-spinning and spinning black holes (for generic 
orbits and arbitrary spin orientations), and that the test- 
particle limit is recovered non-perturbatively, (i.e., at all 
PN orders). This is because in order to build our mod- 
els, in Sec. II C we started from the Hamiltonian derived 
in Ref. [59], which correctly reproduces the Mathisson- 
Papapetrou-Pirani equation describing the motion of a 
classical spinning particle in a generic curved space- 
time [63-67]. The EOB models that we present in this 
paper share this feature with our earlier model [46] , which 
was valid through 3PN order in the non-spinning sec- 
tor and through 2.5PN order in the spinning sector, but 
not with other EOB models for spinning black-hole bi- 
naries, which recover the test-particle limit only approx- 
imately [15]. 

Other attractive features of our models are evident 
when considering configurations with spins parallel to 
the orbital angular momentum, which correspond, in the 
effective EOB dynamics, to a particle moving on equa- 
torial orbits. For aligned spins and equatorial orbits, in 
fact, both the models with dynamical and non-dynamical 
spin mapping predict the existence of an innermost sta- 
ble circular orbits (ISCO), for all values of the system's 
parameters. This feature is again shared by our earlier 
model [46], but not by other EOB models for spinning 
black- hole binaries [15], which do not present ISCOs for 
large values of the spins. While the non-existence of an 
ISCO is not necessarily a sign that a model is flawed, 
its presence helps reproduce the results of numerical- 
relativity simulations for binaries with aligned spins [68] . 

To calculate the radius and the orbital angular mo- 
mentum at the ISCO for our EOB models, we solve nu- 




FIG. 1: The spin parameter of the binary at the ISCO given 
by Eq. (67) for the 2.5PN and 3.5PN EOB models with dy- 
namical mapping of the spins, for binaries having spins par- 
allel to L, mass ratio q = 7712/mi and spin-parameter projec- 
tions onto the direction of L given by \i = X2 — X- 



merically the following system of equations 

(r,p r = 0,L Z ) 



o rj-improved ( 



dr 
d 2 HZr vcd (r,p r =0,L z ) 



real 



dr 2 



= 0, 
= 0, 



(64) 
(65) 



with respect to r and L z = p$. The solutions can then 
be used to evaluate the ISCO frequency via 



o 



f)ff improved/ _ n j isco n 



dL z 



(66) 



which follows immediately from the Hamilton equations. 

The values of r lsco and L z can also be used to 

calculate the binding energy at the ISCO via -Ebind = 

^roa? r ° VC — M ■ This quantity is interesting because it 
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FIG. 2: The same as in Fig. 1 but for the binding energy of 
the binary at the ISCO. 



FIG. 3: The same as in Fig. 1 but for the ISCO frequency. 



corresponds to the mass lost in gravitational waves dur- 
ing the binary's inspiral, and is therefore a lower limit to 
the total mass loss, to which it reduces for r\ — > (when 
the fluxes during the merger and the ringdown become 
negligible [46]). Similarly, one can estimate the spin of 
the binary at the ISCO via 



03 



Xik 



Sf 



S 2 + ^sco 



(M + E^y 



(67) 



This expression clearly neglects the mass and angular 
momentum lost during the merger and ringdown phases, 
but it is useful as qualitative diagnostics of our model, 
and it reduces to the spin of the final black-hole remnant 
when r\ — > (again, because in this limit the fluxes during 
the merger and the ringdown become negligible [46]). 

We re-write the metric potentials A t and A r given in 
Eqs. (33), (34), using the "log-model" of Ref. [46] [see 
Eqs. (5.71) and (5.73)-(5.83) of that paper], and assume 



K{rj) = 1.447 -0.1574 77- 9.082 t? 2 



(68) 



The value of K(rf) for r\ = ensures [46, 69] that the 
ISCO frequency for extreme mass-ratio non-spinning bi- 
naries predicted by our EOB models agrees with the ex- 
act result of Rcf. [70], which calculated the shift of the 
ISCO frequency due to the conservative part of the self- 
force. The linear and quadratic terms in rj in Eq. (68) arc 
such that our EOB models accurately reproduce numer- 
ical relativity simulations for non-spinning binaries with 
mass ratios ranging from q = 1/6 to q = 1 [68]. 

We fix the gauge parameters to the following values: 



a o = -2 7?! 



bo = -p, 



(69) 



Ol 



« 2 



6i 



16 



V (9 + 5r?) . 



77(7-877), & 2 = g/? (17 - 577) 



16 



V 



:T 



(70) 



which wc determine by requiring that all the terms in- 
volving A r p-n/S cancel out in A^» and A^, [Eqs. (51), 
(52)], or equivalently in Qs2.5PN and Q s 3. 5PN [Eqs. (62), 
(63)] . Different choices of the gauge parameters produce 
qualitatively similar results for the ISCO quantities that 
we described above. 

Focusing on systems with spins aligned with the or- 
bital angular momentum _L, and denoting with Si. 2 = 
Xi,2 ™? 2 the projections of the spins along the direc- 
tion of L, we consider binaries with Xi = X2 = X an d 
mass ratios q = 7772/7711 = 0.1, 0.5 and 1. In particular, 
in Figs. 1-3 we show how the ISCO quantities described 
above change as a consequence of including the 3.5PN SO 
terms in our EOB model with dynamical spin mapping. 
More specifically, we calculate Ci lsco M, E^^/M and 
Xisco using the Hamiltonian (54), with and without the 

(2) 
3.5PN terms given by A^, . As can be seen, the inclusion 

of the 3.5 PN terms does not change the ISCO quantities 

significantly for \ < 0, while small differences appear for 

X > 0. (In the case of f2 ISCO M, however, these differences 

grow quite large when \ — > 1.) Overall, Figs. 1-3 suggest 

that the model has reasonable convergence properties for 

radii r > r lsco . 

The results for the model with non-dynamical spin 
mapping are similar [i.e., a comparison of the ISCO quan- 
tities calculated using the Hamiltonian (58), with and 
without the 3.5PN term Q3.5PN, gives similar results]. 
In general, however, the model with non-dynamical spin 
mapping presents lower values for J7 ISCO M at high spins 
and for comparable mass ratios (see Fig. 4, where we 
compare the 3.5PN models with dynamical and non- 
dynamical spin mapping). 

Another attractive feature of our models is the exis- 
tence of a peak of the orbital frequency during the plunge 
starting at the ISCO. More precisely, we assume that the 
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FIG. 4: The ISCO frequency for the 3.5PN EOB models with 
dynamical (dyn) and non-dynamical (non-dyn) mapping of 
the spins, for binaries having spins parallel to L, mass ratio 
q = rri2/mi and spin-parameter projections onto the direction 
of L given by \i = X2 = X- 




FIG. 5: The same as in Fig. 1, but for the maximum of the 
orbital frequency during the plunge. 



effective particle starts off with no radial velocity at the 
ISCO (thus having angular momentum L lsco and energy 
E ISOO ), and we evolve the geodesic equations by calculat- 



ing the radial momentum p r during the plunge from en- 
ergy and angular momentum conservation. We then cal- 
culate the orbital frequency £1 = dH™? rove /dL z along 
the trajectory and find that it presents a peak r2 max . 
This is not surprising because the same behavior was 
observed to be generic in our earlier model [46]. The 
values of MQ max for binaries with spins parallel to L, 
as function of x — Xi = Xi-, are shown in Fig. 5 for 
mass ratios q = 1,0.5 and 0.1, for the EOB model with 
dynamical spin mapping at 2.5PN and 3.5PN. As can 



FIG. 6: The same as in Fig. 4, but for the maximum frequency 
during the plunge. 



be seen the differences introduced by the 3.5 PN terms, 
although reasonable, are larger than for the ISCO quan- 
tities. This may be because the plunge happens at radii 



that are smaller than 



and approach the horizon's 



radius, thus making the higher order PN terms more and 
more important. The results for the model with non- 
dynamical spin mapping are generally similar, although 
they differ slightly at high spins. In particular, in Fig. 6 
we compare the 3.5PN models with dynamical and non- 
dynamical spin mapping. As can be seen, for q = 0.5 and 
q = 1 the predictions of the two models are very close, 
while for q — 0.1 the model with dynamical spin mapping 
presents somewhat lower maximum frequencies. 

Also, we stress that the values of Mf2 max for spin an- 
tialigned with the angular momentum (i.e., \i = X2 = 
X < 0) are quite sensitive to the values of the gauge pa- 
rameters do— as and 6o _ ^3- For instance, setting all the 
gauge parameters to makes the behavior of MO max 
with x non-monotonic if the 3.5PN models (both with 
dynamical and non-dynamical spin mapping) are consid- 
ered. This effect does not appear in the 2.5PN models, 
and can in principle be important for the calibration of 
our model with numerical-relativity simulations. More 
details on this will be given in a follow-up paper [68]. 
Even worse, when the gauge parameters are set to zero 
the difference in Mf2 max between the 2.5PN and 3.5PN 
models is larger than in Fig. 5, a sign that the model 
probably converges more slowly in this gauge. In light 
of this, it seems preferable to use the gauge parameters 
(69)-(70), which by canceling out the radial momentum 
A r p ■ n/T, from A^V and A^.« (and from Qs2.5PN and 
Ss3.5Pn) provide a rather regular and monotonic behav- 
ior for MO max and reasonable differences between the 2.5 
and 3.5PN models. 

Finally, in Fig. 7 we show the predictions of our EOB 
model with dynamical spin mapping for the ISCO fre- 
quency of a system with q = ms/mi = 10~ 3 , xi — X 
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FIG. 7: The shift of the ISCO frequency cn, defined in 
Eq. (71), for the 2.5PN and 3.5PN EOB models with dynam- 
ical mapping of the spins, for a binary having spins parallel 
to L, mass ratio q = 7712 /mi = 10 - and spin-parameter pro- 
jections onto the direction of L given by xi = X an< i X2 = 0. 



and \2 = (the results for the EOB model with non- 
dynamical spin mapping are similar) . More precisely, we 
show the fractional deviation from the Kerr ISCO fre- 
quency normalized by the mass ratio, 



1 ( ^.sco^l* 

cn = - 



q \n 



30-^lKerr 



1 



(71) 



as a function of x, as proposed in Ref. [69]. This ISCO 
shift is caused by the conservative part of the self-force 
and has been calculated exactly by Ref. [70] in the case 
of a Schwarzschild spacetime (\ = 0). The results of 
Ref. [70] is c n = 1.2513 + 0(g) [see also Ref. [71]], and is 
denoted by a filled circle in Fig. 7. As can be seen, both 
the 2.5 and 3.5PN models predict cq > 0, except when 
X <) 0.83. This change of behavior of the EOB predic- 
tion is common also to our earlier model of Ref. [46] , and 
might have important implication for configurations that 
might violate the Cosmic Censorship Conjecture [72, 73]. 
However, the behavior of cn, which seems to diverge as x 
approaches 1, suggests that this might simply be a spuri- 
ous effect due to the incomplete knowledge of the function 
K [Eq. (68)] and to the fact that the EOB model only re- 
produces the SS coupling at leading PN order (2PN). As 
mentioned in Ref. [46] , K may in general depend not only 
on 77 but also on % 2 , and these spin-dependent terms can 
be very important for near-extremal spins, and so will 
the 3PN SS couplings. 

It is therefore possible that after reconstructing the full 
functional form of K (by comparing to future self-force 
calculations in Kerr or to numerical-relativity simulations 
for spinning binaries) and extending the EOB model to 
include the 3PN SS couplings, cq might remain positive 
even at high spins. 



IV. CONCLUSIONS 

Recently, Ref. [19] has computed the 3.5PN SO effects 
in the ADM Hamiltonian. We have taken advantage of 
this result and extended the EOB Hamiltonian of spin- 
ning black holes to include these higher-order SO cou- 
plings. 

Building on previous work [39, 58], and in particular 
on the EOB Hamiltonian of Refs. [46, 59], which repro- 
duces the SO test-particle couplings exactly at all PN 
orders, we have worked out two classes of EOB Hamil- 
tonians, which differ by the way the spin variables are 
mapped between the effective and real descriptions. One 
class of EOB Hamiltonians is the straightforward exten- 
sion to the next PN order of the EOB Hamiltonian of 
Ref. [46]. It uses a mapping between the real and effec- 
tive spin variables that depends on the dynamical orbital 
variables p 2 , n • p and r. By contrast, the other class of 
EOB Hamiltonians uses a mapping between the real and 
effective spin variables that does not depend on these 
dynamical orbital variables. We achieved this result at 
the cost of modifying the Hamilton- Jacobi equation of a 
spinning test-particle. 

Quite interestingly, when restricting to spins aligned 
or antialigned with the orbital angular-momentum and 
to equatorial circular orbits, we find that the predictions 
of these two classes of EOB Hamiltonians for the ISCO 
frequency, energy and angular momentum, and for the 
maximum of the orbital frequency during the plunge are 
generally similar. However, for high spins the model with 
dynamical mapping of the spins may present somewhat 
lower maximum frequencies and larger ISCO frequencies. 

As pointed out originally in Ref. [58], several gauge 
parameters can enter the canonical transformation that 
maps the real and effective Hamiltonians. If the Hamil- 
tonian were known exactly, i.e., at all PN orders, then 
physical effects should not depend on these parameters. 
However, since we know the Hamiltonian only at a cer- 
tain PN order, we expect these gauge parameters to lead 
to non-negligible differences. In fact, we obtained that 
when setting all the gauge parameters to zero, the max- 
imum frequency during the plunge has a non-monotonic 
dependence on the spins, and varies quite significantly 
as a consequence of the inclusion of the 3.5 PN SO cou- 
plings. We found instead that when choosing the gauge 
parameters so that the terms depending on the radial mo- 
mentum disappear from our spin mapping (in the model 
with dynamical spin mapping) or from the modifications 
to the Hamilton- Jacobi equation (in the model with non- 
dynamical spin mapping), the maximum frequency dur- 
ing the plunge has a much more regular behavior and 
varies by small amounts when adding the 3.5PN SO cou- 
plings. This suggests that such a choice of the gauge pa- 
rameters may accelerate the convergence of the model's 
results in the strong-field region where the plunge takes 
place. 

The EOB Hamiltonians derived in this paper can be 
calibrated to numerical-relativity simulations with the 
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goal of building analytical templates for LIGO and Virgo 
searches. A first example was obtained in Ref. [47], where 
the EOB Hamiltonian at 2.5PN order in the SO couplings 
of Ref. [58] was calibrated to two highly-accurate numer- 
ical simulations. Results that use the EOB Hamiltonian 
at 3.5PN order developed in this paper will be reported 
in the near future [68]. 

Lastly, while finalizing this work, Ref. [74] appeared 
in the archives as a preprint. Both this paper and 
Ref. [74] derive the effective gyromagnetic coefficients 
[see Eq. (31)], but with two different methods. Our 
computation uses the Lie method to generate both the 
purely-orbital and the spin-dependent canonical transfor- 
mations, while Ref. [74] first applies explicitly the purely- 
orbital transformation from ADM to EOB coordinates, 
and then uses Eq. (7) to account for the effect of a spin- 
dependent canonical transformation. As a result of these 
different procedures, and as discussed in Sec. II A, the 
2.5PN gauge parameters in our spin-dependent canoni- 
cal transformation coincide with those of Ref. [74], but 



the 3.5PN gauge parameters have different meanings in 
the two approaches and therefore do not coincide. How- 
ever, by suitably expressing our 3.5PN gauge parameters 
in terms of those of Ref. [74], we find that our effec- 
tive gyromagnetic coefficients fully agree with those of 
Ref. [74] . This amounts to saying that our gyromagnetic 
coefficients agree with those of Ref. [74] up to a canonical 
transformation, and are therefore physically equivalent. 
More importantly, in this paper we have focused on and 
worked out two classes of EOB Hamiltonians that are 
different from the one considered in Ref. [74]. 
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